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ABSTRACT 
The overuse of B-lactam antibiotics has caused drug-resistant bacteria, including 
NDM-1, a Metallo-B-lactamase that renders inhibitors ineffective. NDM-1 is 
in Klebsiella 
Developing an NDM-1 inhibitor is vital for managing bacterial resistance and 


found pneumoniae, causing worldwide resistant infections. 
public health. This study investigated L-aspartic acid B-benzyl ester 7-amido-4- 
methylcoumarin (Asp (OBzl)-AMC) as a potential NDM-1 inhibitor through in 
vitro and in silico evaluations. Enzyme kinetics indicated competitive inhibition, 
with an ICs0 value of 30.4 + 5.0 uM and a Ki value of 11.4 + 1.2 uM. Molecular 
docking showed Asp (OBzl)-AMC forming hydrogen bonds with NDM-1's active 
site residues. Asp (OBzl)-AMC displayed favorable ADME properties, making it 
a promising drug candidate. The study demonstrates that Asp (OBzl)-AMC 
effectively binds to NDM-1 and inhibits its activity. The compound's attachment 
to the naphthalene ring is crucial for protecting B-lactam antibiotics from NDM-1 
damage. These findings suggest that Asp (OBzl)-AMC holds potential as an 


NDM-1 inhibitor to combat bacterial resistance and maintain public health. 


Keywords: Metallo-beta-lactamase, enzyme, inhibition, kinetics, molecular 


dynamics 


1. INTRODUCTION 


Inhibition of bacterial infections is essential for maintaining human health (Al- 
Jalali and Zeitlinger, 2018; Shamriz and Shoenfeld, 2018). 6-lactam antibiotics, 
including penicillins, carbapenems, and cephalosporins, are widely employed to 
treat such infections (Edwards, 2000; Essack, 2001; Liu et al., 2018; Paterson et al., 
2001). However, excessive use of these antibiotics has led to the emergence of 
drug-resistant bacteria, resulting in a decline in the efficacy of these drugs and 
raising concern for public health (Khan-Ali et al., 2017; Khan-Maryam et al., 2017). 
The primary cause of antibiotic resistance is the proliferation of f-lactamases, 
enzymes that catalyze the hydrolysis of B-lactam antibiotics. These enzymes are 
classified into four classes, A through D, with class B requiring zinc ions in the 


1 of 10 


ANALYSIS ARTICLE OPEN ACCESS 


active site for activity, and are referred to as metallo-B-lactamases (MBLs) (Hawk et al., 2009). 

The subclass B of Metallo-B-lactamases (MBLs) is further classified into three subcategories: B1, B2, and B3. Currently, f- 
lactamase inhibitors, including clavulanic acid, sulbactam, and tazobactam, are commercially available and can be used in 
conjunction with antibiotics to treat some antibiotic-resistant infections, but all are ineffective against MBLs. Despite the presence of 
potential inhibitors for MBLs (Crowder et al., 2006; Materon et al., 2003), none are commercially available for use against these 
enzymes. Given the prolonged drug development process, it is imperative to identify effective inhibitors for MBLs. 

A strain of Klebsiella pneumoniae emerging from India has been found to produce an MBL of the B1 subclass, known as New 
Delhi Metallo-B-lactamase-1 (NDM-1) (Brem et al., 2016; Groundwater et al., 2016). This enzyme contains two zinc ions, Zn1 and 
Zn2, which are directly involved in its catalytic function. The presence of NDM-1 has since been reported in various regions 
worldwide (Heinz et al., 2019; Johnson and Woodford, 2013; Tran et al., 2015). NDM-1's resistance to 6-lactamase inhibitors is 
primarily due to its ability to hydrolyze the amide bond of the B-lactam ring (King et al., 2012). The resistance of pathogenic bacteria 
to treatment presents a significant challenge to disease management (Ejaz et al., 2020). As such, targeting NDM-1 may potentially 
inhibit the resistance of pathogenic bacteria. 

Genes associated with NDM-1 can be transferred to Escherichia coli (Shah et al., 2015; Tuem et al., 2018; Zhang et al., 2017). In 
recent years, resistance to E. coli has increased among the population (Wang et al., 2015). NDM-1 plays a significant role in those 
bacterial resistance mechanisms (Du et al., 2017). Plasmids containing NDM-1 in patients can spread among different species of 
Enterobacteriaceae, another common pathogenic bacteria in healthcare settings (Martino et al., 2019). Bacterial infections resulting 
in this type of resistance to antibiotics not only increase medical expenses but also poses a threat to human health (Yang et al., 2017). 

Therefore, it is of paramount importance to develop an effective drug that can inhibit NDM-1. In this study, we aimed to 
investigate the potential of L-aspartic acid 6-benzyl ester 7-amido-4-methylcoumarin (hereafter referred to as Asp (OBzl)-AMC) as 
an NDM-1 inhibitor, given its possession of a potential zinc-binding group, ester and amide. The ability of the compound to 
effectively bind the zinc ion in the active site was evaluated using both in vitro and in silico methods. Enzyme kinetics inhibition 
studies and a series of molecular dynamics simulations of NDM-1 in complex with the compound Asp (OBzl)-AMC were 
conducted. These results may offer a promising approach for developing a potential drug that targets NDM-1. 


2. MATERIALS AND METHODS 


General Procedures 

All chemicals were obtained through Sigma or other quality manufacturers. The NDM-1 gene from K. pneumoniae was inserted into 
the pET100 vector using the Champion pET Directional TOPO Expression Kit from Invitrogen (Waltham, MA, USA). The enzyme 
was purified using previously established methods with the pET100-NDM-1 construct (Schlesinger et al., 2011). The product's 
purity was determined through Ni® chromatography using 10% SDS-PAGE. Additional purification was achieved through gel- 
filtration chromatography if necessary. The enzyme was stored at -20°C in 50 mM 3-(N-morpholino) propane sulfonic acid (MOPS) 
(pH 7.0), 0.050mM ZnSO4, and 30% (v/v) glycerol for further experimentation when it was determined to be 95% pure or higher. 
Asp (OBzl)-AMC was used as a substrate and was prepared by dissolving it in dimethyl sulfoxide (DMSO) before being diluted to 


various concentrations for enzyme activity assays. 


Inhibition Tests 

The purified NDM-1 enzyme activity was assessed by adding it to solutions of nitrocefin (ranging from 10 to 70 UM) in 50 mM 
MOPS (pH 7.0) at a final concentration of 0.13 g/mL, in a total volume of 1 mL, using quartz cuvettes at room temperature. The 
increase in absorbance at 485 nm resulting from the hydrolysis of nitrocefin (¢ = 15,900 M-cm"!) was used as a measure of enzyme 
activity, and the assays were conducted in triplicate (Fast et al., 2001). The ICs0, or the concentration of inhibitor required to inhibit 
50% of enzymatic activity, was determined using an enzyme assay with Asp (OBzl)-AMC as the inhibitor and nitrocefin as the 
substrate. The assay was performed using a fixed concentration of nitrocefin (0.4 mM) and a fixed enzyme concentration (0.13 
ug/mL or 4.8 nM) after a 1-minute incubation period with varying inhibitor concentrations. To further analyze the mode of 
inhibition, a separate series of assays were conducted at fixed inhibitor concentrations (0 — 0.030 mM) with varying substrate 
concentrations (0.1 to 0.4 mM) and analyzed using Lineweaver-Burk plots and non-linear regression through Sigma Plot version 
11.0. The competitive inhibition equation, v = Vmax:S/[Km-(1 + 1/Ki) + S], was used for data analysis. All experiments were performed 


in triplicate. 
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Molecular Docking Simulations 

This study used the crystal structure of the Metallo-B-lactamase NDM-1 (PDB id: 4EY2) as a starting point. The structure of Asp 
(OBzl)-AMC was computationally generated using the Avogadro program, with energy minimization carried out using the General 
Amber Force Fields (Hanwell et al., 2012). A molecular docking simulation was performed for each ligand using the Autodock Vina 
program, with full flexibility for the ligand and a rigid receptor. The active site of NDM-1 was defined as the center of a grid box 
with dimensions of 45 A x 45 A x 45 A, with the x-y-z coordinates of (-0.509, 8.133, 23.576). The lowest energy docking results for 
each ligand were used for further analysis through molecular dynamics. The Protein-Ligand Interaction Profiler (PLIP) online 
server was utilized to analyze the interactions from the Autodock results and the final conformations obtained from molecular 


dynamics (Adasme et al., 2021). 


Molecular Dynamics Simulations 

All molecular dynamics simulations were conducted using the Gromacs program version 2021.4. The AMBER99SB-ILDN force 
field, modified with recently published zinc (II)-binding residues parameters, was employed in these simulations (Macchiagodena 
et al., 2019; Macchiagodena et al., 2020). The modifications to the force field were implemented to ensure accurate coordination of 
the active site zinc ions, thereby enhancing the stability of the metalloprotein during the dynamics simulation. The SPC/E water 
model was used in the simulation. To create the topologies of the ligands for use in Gromacs, the Acpype program was utilized 
(Sousa-da-Silva and Vranken, 2012). The initial step of each dynamics simulation involved solvating the system and neutralizing it 
by adding three chlorine ions. An energy minimization was then performed for 50,000 steps using the steepest descent algorithm. 
The system was equilibrated in two phases: A 100 ps run using the NVT ensemble with the leap-frog integrator, which stabilized 
the temperature at 300 K, and a 100 ps run using the NPT ensemble with the leap-frog integrator, which stabilized the pressure at 1 
bar. 

Position restraints were applied to the receptor and ligand during each equilibration run. The Berendsen V-rescale thermostat 
was used to regulate the temperature, and the Parrinello-Rahman barostat was used to regulate the pressure. The final production 
MD simulation was run for 100 ns, with the position restraints removed. Two replicate simulations were conducted for the NDM-1: 
Asp (OBzl)-AMC complex system with different initial velocities and were compared to identify consistent behavior. The replicate 
with the best fit was then used to further analyze and compare inhibitors. The root-mean-square deviation (RMSD) of each inhibitor 
trajectory relative to the protein backbone was calculated using the RMSD module of Gromacs. Additionally, the protein backbone 
RMSD relative to the energy-minimized conformation was calculated for each system to provide insight into the protein stability 
throughout the simulations. 


Free Energy Calculations 

Each simulation's average end-state free energy was calculated using the molecular mechanics with the generalized Born and 
surface area (MM-GBSA) model implemented in the gmx_MMPBSA program (Valdés-Tresanco et al., 2021). To perform this 
calculation, 100 snapshots were taken from each simulation. The binding free energy and a per-residue decomposition of the free 
energy were calculated according to a previously published method (Huckleby et al., 2022). 


3, RESULTS 


ICso value determination 

To investigate the potential inhibitory effect of Asp (OBzl)-AMC on NDM-1 activity, we conducted an assay to determine its ICso 
values. Various Asp (OBzl)-AMC concentrations, ranging from 0.5 UM to 100 mM, were used to obtain the ICs0 values. The obtained 
data were analyzed using the concentration-response equation vi/Vo = 1/(1 + ([I]/ICs0)'), where I denote the inhibitor and h is the Hill 
coefficient (values of h were between 0.5 and 1). The data points of Asp (OBzl)-AMC exhibited a satisfactory fit to a semilogarithmic 
concentration-response plot, as demonstrated in Figure 1. By calculating the ICs0 value at 50% inhibition, we determined that the 
ICso of Asp (OBzl)-AMC was 30.4 + 5.0 uM. 
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Figure 1 Concentration-response plots for the inhibition of NDM-1 by Asp (OBzl)-AMC. The letter I denote the inhibitor Asp 


(OBzl)-AMC. The data points are based on the average values obtained from triplicate analysis. 


Mode of inhibition assay 

To investigate the mode of inhibition, inhibitory enzyme assays were performed using varying concentrations of the substrate 
cephalexin and Asp (OBzl)-AMC. The Lineweaver-Burk plots for Asp (OBzl)-AMC are presented in Figure 2, demonstrating an 
intersection of lines at the y-axis. The Ki values were calculated from the slopes of the plots, resulting in a Ki value of 11.4 + 1.2 uM. 


The apparent Km values were determined from Figure 2 to calculate the Ki values. 
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Figure 2 Lineweaver-Burk plots of the inhibition of NDM-1 by Asp (OBzl)-AMC. Cephalexin was used as the substrate. The assays 
were carried out in a buffer of 50 mM MOPS at pH 7.0. Kinetic constants were determined by fitting the data to the equation v = 
Vimax'S/(Km-(1 + 1/Ki) + S), which indicates competitive inhibition. The concentrations of the inhibitor Asp (OBzl)-AMC used were 0 
(orange), 20 (blue), and 40 uM (gray). All data points are based on the average values obtained from triplicate analysis. 


Molecular docking simulations 
Molecular docking simulations were employed to simulate the possible inhibition of NDM-1 by Asp (OBzl)-AMC. The complex 
with the lowest energy and a AG value of -7.8 kcal/mol was selected for further dynamics analysis. The conformation of the chosen 
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complex demonstrated that the carboxyl group interacts with a Zn (II) ion, while a potential hydrogen bond was observed between 
the compound and active site residues His249, as illustrated in Figure 3. 


HIS249 


Figure 3 Autodock binding conformation selected for dynamics. The compound Asp (OBzl)-AMC was centered in pink color. The 
oxygen element was colored in red and other elements were colored in purple or dark purple. 


We employed a PLIP to comprehensively analyze the protein-inhibitor interaction. The results, depicted in Figure 4, reveal the 
occurrence of hydrophobic, hydrogen bonding, and salt bridge interactions between the inhibitor and NDM-1. A detailed account 
of the chemical interactions can be found in Table 1. 


TRP92 


ASN219 


Figure 4 Evaluating docking results with PLIP. The protein is depicted in red, the inhibitor in orange, and the charge center in 
yellow. Hydrophobic interactions are shown by gray dashed lines and hydrogen bonds by purple solid lines. 


A 100 ns molecular dynamics simulation was conducted to investigate the molecular interactions between the ligand and 
receptor in the context of hydration. The lowest energy docking result was used to initiate the simulation, which resulted in a AG 
value of -7.8 kcal/mol. Although the initial molecular docking suggested a hydrogen bond between His249 and the ligand, the 
molecular dynamics simulation showed that Asn192 and His222 also formed hydrogen bonds due to equilibrium processes. The 
simulation snapshots captured in Figure 5 revealed the hydrogen interactions at various stages of the simulation, including between 
the naphthalene ring of the inhibitor and the epsilon-amino group of Lys183, as well as between the naphthalene ring of the 
inhibitor and the carbonyl oxygen of the peptide of Leu190. These observations suggest that the naphthalene ring moiety of the 
inhibitor is critical for its interaction with the target enzyme. 


Drug Discovery 17, e28dd1945 (2023) 5 of 10 


ANALYSIS ARTICLE OPEN ACCESS 


Table 1 The non-covalent interactions (hydrogen bonds and hydrophobic interactions) between the inhibitor and the enzyme. 


Hydrogen Bonds 

feds Ronen: AA Distance | Distance | Donor Protein cath 
H-A D-A Angle Donor Chain 

il 219 Asn 2.06 3.06 173.57 v v 

2 249 His 2.26 2.86 116.50 x v 

Hydrophobic Interactions 

Index Residue AA Distance 

1 34 Tle 3.52 

2 34 Tle 3.50 

3 69 Phe 3.22 

4 72 Val 3.83 

5 72 Val 3.98 

6 92 Trp 3.79 

7 219 Asn 3.90 


The RMSD values for the atomic positions obtained from the simulation were analyzed to examine conformational binding 
variations in a hydrated environment (Figure 6). The RMSD of the peptide backbone of the enzyme remained stable after 0.7 ns, 
indicating that the molecular system was well-behaved throughout the simulation. Similarly, the RMSD of the inhibitor atoms 
stabilized after 0.2 ns, indicating that the system was well-equilibrated. 


Md = mh . ual 
Figure 5 Snapshots of the molecular dynamics simulation for Asp (OBzl)-AMC with NDM-1 at the beginning of the simulation (A) 
and 75 ns (B). 


To gain a more detailed understanding of the interactions between NDM-1 and Asp (OBzl)-AMC, the decomposition data for 
the system were analyzed (Figure 7). The MM-GBSA technique was employed to determine the binding free energy of Asp (OBzl)- 
AMC in complex with NDM-1, with a AG binding value of -21.78 kcal/mol estimated based on the inhibitor confirmation obtained 


from the molecular dynamics simulation. 


Pharmacokinetics properties 

The pharmacokinetics properties of a lead compound, commonly referred to as Absorption, Distribution, Metabolism, and 
Excretion (ADME), are crucial for a drug's success. These parameters track the drug's behavior in the body from administration to 
excretion through sweat, urine, or feces (Dowty et al., 2010). The distribution volume to tissues and target sites affects the drug's 
bioavailability and minimizes side and toxic effects. Evaluating a drug candidate's efficacy and toxicity before clinical trials can 
reduce the risk of unwanted effects in human or animal models. The online Swiss ADME server can assist in reducing costs by 
predicting potential drug rejection. Our study evaluated the ADME properties, including lipophilicity, water-solubility, drug- 
likeness, and medicinal chemistry, of Asp (OBzl)-AMC. The ADME properties of Asp (OBzl)-AMC were found to be in an optimal 
range, as listed in Table 2. 
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Figure 6 The RMSDs for the complexes and the compounds. (A) NDM-1 backbone RMSD from the NDM-1: Asp (OBzl)-AMC 
complex (B) RMSD of Asp (OBzl)-AMC from the NDM-1: Asp (OBzl)-AMC complex. Plotted as running averages for visual clarity. 
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Figure 7 Per-residue energy decomposition diagrams. (A) Binding free energy decomposition of the NDM-1: Asp (OBzl)-AMC 


complex, (B) Per-residue energy decomposition heatmap. 


Table 2 Identification of Physical, Chemical, Pharmacokinetic, and Drug-Likeness Properties of Asp (OBzl)-AMC. 


Drug Discovery 17, e28dd1945 (2023) 


Properties Parameter Asp (OBzl)-AMC 
MW (g/mol) 290.27 g/mol 
Num. heavy atoms 21 
Num. arom. heavy atoms 10 
Physicochemical properties Num. rotatable bonds 5 
Num. hydrogen bond acceptors | 6 
Num. hydrogen bond donors 3 
Molar refractivity 75.86 
Lipophilicity Log Pow 0.58 
Water solubility Log S (ESOL) -0.61 
GI absorption High 
Pharmacokinetics Blood-Brain Barrier permeant No 
CYP3A4 inhibitor No 
: Lipinski violation 0 violation 
Drug likeness : ee 
Bioavailability score 0.55 
Medi. chemistry Synthetic accessibility 3.17 
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4. DISCUSSION 


We have demonstrated the potential for inhibition of NDM-1 by an L-aspartic acid derivative that can bind to the zinc ion in the 
active site. The effectiveness of the zinc ion binding was investigated by way of in silico and in vitro analyses. The MD simulations of 
NDM-1 in complex with the L-aspartic acid derivative revealed the details of the binding interactions. In addition, kinetic analyses 
were performed to examine the inhibitory potential of each compound. The performance of the MD simulations for the NDM-1: 
Asp (OBzl)-AMC complex clearly showed the ability of the L-aspartic acid derivative to remain bound to the zinc in the active site. 
The naphthalene ring moiety of Asp (OBzl)-AMC appears to confer unique advantages in interacting with the Lys183 and Leu190 
side chains in various ways. Also, according to PLIP and per-residue decomposition data obtained from MD simulations, the 
hydrophobic interactions around the aromatic ends of Asp (OBzl)-AMC contributed significantly to stabilizing the compound 
within the NDM-1 active site. 

Furthermore, it should be mentioned here that the molecular docking results appeared to be adjusted for the beginning of the 
molecular dynamics simulation, as the target enzyme is rigid in Auto Dock, while the latter offers total flexibility for the complex. 
The experimental and in silico inhibition studies suggest that the inhibition mode is likely to be competitive. These compounds bind 
strongly to the active site zinc ions, as evidenced by the MD simulations, reaffirming their potential as competitive inhibitors. 
Moreover, the ICs0 and Ki values were in the low micromolar range. Given all these findings, Asp (OBzl)-AMC could serve as a 
promising drug candidate and a lead compound for further development. Structural modifications can be made to enhance the 
drug's effectiveness, and a structure-activity relationship analysis should be conducted. It is also crucial to evaluate the cytotoxicity 
and specificity of these compounds before using them in clinical settings. 


5. CONCLUSION 


To summarize, we have created an L-aspartic acid derivative with a strong affinity for NDM-1, a major contributor to B-lactam 
antibiotic resistance. The compound has demonstrated potential for inhibiting NDM-1, with ICs0 values in the micromolar range 
and a Ki value in the low micromolar range indicating strong binding. Molecular dynamics simulations have provided insight into 
the specific amino acids involved in the complex formation between Asp (OBzl)-AMC and NDM-1. The experimental inhibition 
studies have supported the idea that the mode of inhibition is competitive, which aligns with data obtained from in silico analysis. 
The attachment of the L-aspartic acid derivative to the naphthalene ring is believed to be crucial for interaction with zinc ions, 
contributing to the compound's effectiveness in inhibiting NDM-1. We anticipate that this compound will protect B-lactam 
antibiotics from the damage caused by NDM-1. 
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